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The dynamics of a pair of coupled harmonic oscillators in separate or common 
thermal environments is studied, focusing on different indicators of quantumness, such as 
entanglement, twin oscillators correlations and quantum discord. We compare their decay 
under the effect of dissipation and showr, through a phase diagram, that entanglement 
is more likely to survive asymptotically than twin oscillators correlations. 
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1. Introduction 



The characterization of correlations of a quantum state is object of an intense field 
of investigation, due to both its fundamental scientific interest and its importance 
towards the implementation of quantum technologies-L^J Entanglement has been tra- 
ditionally considered as a fundamental resource to obtain quantum computational 
advantages, and has been used as the main indicator of the quantumness of correla- 
tion. Indeed, as shown in Ref. [2J for pure-state computation, exponential speed-up 
only occurs if entanglement grows with the size of the system. Once mixed-state 
computation is considered, however, signatures of quantum speed-up can come out 
using factorized states as, for instance, in the so-called Deterministic Quantum 
Computation with one Qubit (DQCl).'^Decoherence effects due to dissipation are 
known to be detrimental for entanglement that is indeed disappearing after a tran- 
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sient time (and not asymptotically)!^ 

Among many other attempts to quantify qua ntum correlations, a predominant 
role has been assumed by quantum discord!^*^ It has been introduced with the 
aim of capturing all quantum correlations, including entanglement. However, the 
relationship between these two quantities is still unclear, since they seem to capture 
different properties of the states. In Ref. [7] it is shown, for instance, that even if 
discord and entanglement are the same for pure states, mixed states maximizing 
the discord in a given range of classical correlations are actually separable. Recently, 
the analytical expression of quantum discord has been obtained also for Gaussian 
states l^ opening the possibility to use it for continuous variables.!^ 

From a different point of view, the quantumncss of a system can be measured 
through other kinds of indicators, widely developed, for instance, in the field of 
quantum optics. A well-known example are quantu m corr elations between two twin 
beams generated in optical parametric oscillators! I I The quantumness of the 
state of the emitted light is measured by the absence of fluctuations in their in- 
tensities difference. This absence of noise is equivalent to the negativity of this 
variance for normal ordered operators and was first predicted by Reynaud and 
collaborators^ and experimentally measured in Ref. 1131 Our aim in this work is 
to sec if there is any connection between the latter correlations, which we will call 
"twin oscillators correlations" , with entanglement and discord, comparing their de- 
caying and robustness. 

In this paper, we will consider one of the most fundamental interacting sys- 
tems, i.e. two coupled harmonic oscillators, in the presence of dissipation due to 
the interaction with a thermal environment. Two extreme scenarios we are going to 
investigate are represented by the so-called "common bath" , where the two oscilla- 
tors are thought to be so "close" with each other that they interact with the same 
thermal modes, and the case of "separate baths" , where the dephasing channels are 
completely independent.^'^ It was recently show n '^ 1 " I that two identical oscilla- 
tors in the presence of a common bath can exhibit asymptotic entanglement robust 
against decoherence, depending on the bath temperature and initial squeezing. This 
is a very peculiar case and this behavior is generally lost if the two oscillators are 
not identical.LL^ Still, slow decay of entanglement and robust quantum correlations 
appears in the presence of synchronization between detuned oscillators, as shown 
in Ref. [HI 

Studying the dynamics of the system through the master equation approach, we 
want both to analyze the behavior of quantum correlations, considering entangle- 
ment and quantum discord, and classify the global quantumness of the state using 
the variance of the difference of the occupation numbers. First, we will define the 
model and discuss its solution; afterward, the relevant indicators will be defined; as 
a final step, we will study the dynamics of such indicators. 
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2. The model 

Let us consider two quantum harmonic oscihators ahowing for diversity in their 
frequencies and direct couphng as weh as dissipation in thermal environment M^ The 
model describing the system dipped into two identical but independent (separate) 
thermal baths is given by the total H = Hs + Hb + Hsb, where the system 
Hamiltonian 

Hs = Y+2^i^i + Y 2^'^^'^ ^^^'^ (1) 

describes two oscillators with frequencies uJi.2, unitary masses and coupled through 
their positions, 



fe 4=1 \ y 



^--EEPV + ^^r^^ 



is the free Hamiltonian of two (identical) baths of harmonic oscillators (labeled by 
k), and 

HsB^j:>^i'^xi%,+Y^X^^^xi%, (3) 

fc fe 

encompasses the system-bath interaction. 

The case of a common bath is obtained by considering only z = 1 in Hb and 

HcB=^>^kXk{xi+X2). (4) 

fe 

The effective dissipation takes place therefore only in the sum of positions (I'l +X2), 
while in the case of separate baths both positions xi and X2 are independently 
coupled with the thermal bath. 

The thermal bath is assumed hereto be Ohmic, with a Lorcntz-Drude cut-off 
parameter A, and its spectral densitjLLiI is 

J{Q) = -faLu9{A-uj) (5) 

In the following we will consider the cut-off frequency always larger than any fre- 
quency involved in the oscillators free dynamics, that is A >> uji^2 and weak 
coupling 7o. 

The analysis of the dissipation of identical oscillators in common and separate 
baths was given in Refs. [HI [TBI (SQI In these works it is shown that in presence 
of a common bath entanglement can persist in the asymptotic state, in spite of 
dissipation. The master equations describing the time evolution of the reduced 
density matrix of the system (obtained by tracing out the degrees of freedom of the 
bath) when including the frequency diversity of the oscillators have been reported 
in Refs. I17[ 1181 in the weak coupling limit (between system and environment) and 
without relying on the rotating wave approximation. Once the master equation has 
been obtained, we can explicitly write down and solve the equations of motion for 
all the operators moments. If the initial state of the two oscillators is Gaussian with 
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vanishing average positions and momenta, the complete information is contained in 
the matrix of second moments (the covariance matrix). Detailed dynamical equa- 
tions, for the reduced density matrix and for the position and momen ta s econd order 
moments, are given in Ref. [ITl also in the system eigenmodes basia^ and not re- 
produced here. The limit case of high temperatures in presence of a non-Markovian 
environment is studied in Ref. 1211 

3. Quantifiers of quantumness 

In this section, we briefly review the quantifiers we will use to characterize the time 
evolution of our system. 

3.1. Entanglement 

While, in general, measures of entanglement have only been developed for pure 
states, the case of Gaussian density matrices, together with the case of qubits, 
is one of the exceptions, since a necessary and sufficient criterion of separability 
exists.!^ 

For pure bipartite states \4>ab) = X^n "^"l"")!^'")' independently on their na- 
ture, entanglement can be calculated through the von Neumann entropy (en- 
tropy of entanglement) of one of the two reduced density matrices: E = 
-TrA(|(/'AB)(0As|log|(/)As)(0AB|) = -Tr_B (|(/)ab) ('?!)AB | log It/i^B) (0.1S |)- In the 
case of mixed states p, however, the von Neumann entropy cannot be used since 
the mixedness of the reduced density matrices cannot discriminate between entan- 
glement and lack of purity of g. A sufficient criterion (the so-called Peres-Horodecki 
criterion) for detecting entanglement can be obtained by considering the positivity 
of the partial transpose p^^ (or, equivalently, p^^), i.e. of the matri x obta ined by 
only transposing the degrees of freedom of one of the two sub parties ! ^^ Indeed, 
the presence of negative eigenvalues of p^^ witnesses that p has not the form of a 
factorized density matrix. As said before, in the case of Gaussian states, the sep- 
arability of the partial transpose is also necessary to detect entanglement,'^ and 
the modulus of the sum of the negative eigenvalues of p^^ {N) has been shown to 
be an entanglement monotone.'^ 

Since Gausssian states are completely characterized by their first and second 
moments, and first moments can be set to zero with local operations that do no 
modify entanglement, the covariance matrix can be used to check the positivity of 
the partial transpose. The logarithmic negativity, which represents an upper bound 
to the distillable entanglement, is defined as Ej\f = log2(2A/' -I- 1) and is related to 
the smallest symplectic eigenvalue of the covariance matrix of p^^ {^-)'- 

Em = max[0, - loga 2A_]. (6) 

In contrast to other entanglement measures, logarithmic negativity does not reduce 
to entropy of entanglement on pure states. 
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This entanglement quantifier was also eonsidered in the eqntext of coupled dis- 
sipative harmonic oscillators in the mentioned works! ^ ^ < | to | zU | ^l | 



3.2. Quantum discord 

In classical information theory, the mutual information of a bipartite system can 
be calculated through two equivalent formulae related by Bayes rule: we have I{A : 
B) = J{A : B), with I{A : B) = H{A) + H{B) - H{A,B) and J{A : B) ^ 
H{A)~ H{A\B), where H{.) is the Shannon entropy and H{A\B) is the conditional 
Shannon entropy of A given B. 

The quantum counterparts of I(A : B) and J {A : B), however, differ 
substantially!^ By replacing the Shannon entropy with the von Neumann entropy 
of a given bipartite state g [S{g) ~ — Trplogj g], we obtain the quantum mutual 
information 

Iig)^S{gA) + S{gB)~S{g), (7) 

where gA{B) are the reduced states after tracing out party B{A). Due to the nature 
of measures in quantum mechanics, J'{g) depends on the measurement realized on 
B. Classical correlations are then defined as^^ 

Jig){nf} = rniniSigA) - 5(A|{nf })], (8) 

with the conditional entropy defined as S{A\{E^}) = X^jPi5'(g^|£;B), pi = 
TiAsiEf g) and where gAiE^ = Efg/pi is the density matrix after a positive 
operator valued measure (POVM) {{E^}) has been performed on B. Quantum 
discord is defined as the difference between I{g) and J{g)'- 

6A:B{g) = min [S{gB) - S{g) + S{A\{Ef})\ . (9) 

{Bf } 

While the calculation of quantum discord, being based on a minimization procedure, 
is in general an unsolved problem, in the case of Gaussian states, an analytical 
formula has been obtained!^ This allowed to consider quantum discord also in the 
context of continuous variable quantum information!^ 

The quantum discord measures in some sense how much disturbance is caused 
when trying to know about party A when measurin g p arty B, and has been shown 



to be null only for a set of states with measure zerof^^ It was shown to be a useful 
resource in the DQCl algorithm,'^ where the quantum speed up does not rely on 
entanglement and, given its inequivalence to entanglement (except for pure states), 
it hints at more general definitions of what is quantum in a correlation. 

The dynamics of quantum correlations, as quantified by the discord, and mu- 
tual information between quantum harmonic oscillators have been recently studied 
between different oscillators focusing on different parameters regimes and showing 
that the robustness (slow decay) of these correlations is related to the presence of 
a synchronous dynamics!^ 
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3.3. Twin oscillators 

In the context of quantum optics, the discrimination between the predictions of 
classical and quantum theories has been a wide field of investigation. A tool to 
investigate the violation of classical inequalities is given by the variance of the 
difference of the occupation numbers: 

d=(:A2(ni-n2):) (10) 

where, as usual, rii is the occupation number operator of each oscillator and (: 
. :) indicates normal ordering. The quantumness of correlations in the occupation 
numbers derives from the absence of noise when subtracting the oscillators intensity 
fluctuations. This is equivalent (in normal ordering) to the negativity of the variance 
d and is a consequence of the negativity of the Glauber-Sudarshan quasi-probability 
P. 

We note that for identical osciUators Eq. [TOlis identical to (: (ni — 712)^ :) and 
this indicator characterize anti-bunching, being a cross correlation larger than the 
autocorrelation. When the system or the state are not symmetrical in the two com- 
ponents, the negativity of the latter second order moment would be one of many 
possible quantum indicators, implyii ig n egativity of the correspondent P distribu- 
tion not associated to anti-bunching.l^ The existence of these strong correlations 
in optics generally comes from the simultaneous generation of pair s of photons in 
nonlinear processes, and this generally characterizes twin beams.L'^This is a rather 
robust phenomenon in complex spatiotemporal dynamics.'^ 

In the following wc will study the temporal dynamics of d, that we will name 
twin oscillator correlations. Given the Gaussian character of the initial state we w ant 
to study, the fourth order moments can be obtained from the covariance matrix.'^ 
The dynamics of this indicator was already considered in Ref. [17] in comparison 
with entanglement through few examples, showing a similar decay after a finite 
time transient. However, the possibility to get asymptotically twin oscillators was 
not considered there and will be addressed in the next section where we will also 
fully analyze the role of the squeezing of the initial state and of the temperature. 
On the other hand, this indicator (d) has been considered in Rcf . [21] in a different 
regime, for high temperatures focusing on the short time decay in presence of non- 
Mar kovian environment. 

4. Correlations dynamics 

Let us consider an initial two-mode squeezed state 

00 

|*TA/5) = v/T~^E'^"^'l")l")' (11) 

n=0 

where ^ = tanh r and r is the squeezing amplitude. We know that for this state 
d = — 2/i/(l — /^) < 0. We want to study how the different indicators dynamically 
behave considering three different scenarios: (i) the case of different frequencies for 
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a common bath; (ii) the case of equal frequencies for separate baths; (iii) the case 
of equal frequencies (wi = u)2) for a common bath. As shown in Refs. 1151 and I16[ 
for equal frequencies, in the case of a common environment, the asymptotic state 
is expected to be entangled. This is due to the fact that one of the degrees of 
freedom of the system [the mode x^ — {xi — X2)/\^] is actually frozen, since it 
does not interact with the bath, and represents a decoherence-free subspace. Moving 
away from the resonance condition, a full thermalization process takes place, and 
the final (Gibbs) state can be entangled only in the very low temperature regime. 
Motivated by these results, we then want to investigate whether the existence of 
this noiseless channel also protects other aspects of the quantumness of the state 
under evolution. We start A = (no direct coupling between the oscillators) to 
avoid fast oscillations and wi as a scale unit of measure of energy. Then, we choose 
for the system bath coupling 7 = 2 x IO^^/ttwi (weak coupling limit) with cut-off 
A = 20a;i, and squeezing r = 2 in the initial state (fTT|) . 

In Fig. 1, we show the dynamics of logarithmic negativity, quantum discord and 
max[0, — <i/4] {d is scaled by a factor 4 for the sake of comparison). We start consid- 
ering resonant oscillators (wi = CJ2) coupled to separate baths and detuned oscilla- 
tors (wi y^ UJ2) couple d to a common environment. Entanglement and max[0, —d/A] 
vanish in a finite t ime ,Li^ being more fragile than quantum discord, which is expo- 
nentially decayi ng!-^ As a matter of fact, it is known that states with vanishing 
discord are rarej I I and in the presence of dissipation this indicator does not ex- 
perience a sudden death process. 

The case of equal frequencies an d a com mon environment is that one showing 
entanglement in the asymptotic state! I I As we can observe from the right panel 
of Fig. 2, for temperature T = uji, both entanglement and quantum discord reach 
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Fig. 1. Dynamics of logarithmic negativity (gray), quantum discord (black), and max[0, — (i/4] 
(light gray) for T = uj^. Left panel: uj2 = 1.2li;i and common bath; right panel: uJ2 = u)i and 
separate baths. 



July 29, 2011 0:8 WSPC/INSTRUCTION FILE robust 



8 Gian Luca Giorgi, Fernando Galve and Roberta Zambrini 

a stationary regime after a transient phase, while the negativity of d disappears in 
a finite time. By lowering T up to O.lwi, i.e. in a regime where quantum effects are 
stronger, we see that the asymptotic value of entanglement and discord is increased, 
but we also observe that the variance of the difference between the occupation 
numbers becomes negative for infinite time. In Fig. 3, we reproduce the left panel 
of Fig. 2 by adding a finite direct coupling between the two oscillators A = 0.2a;^. 
While the presence of A induces fast oscillations in the three observables, all of them 
are still present in the asymptotic regime. 

Finally, comparing these first figures, it seems then that the sub-poissonian 
character of fluctuations [d < 0) is rather fragile and that oscillators remain twin 
asymptotically only at very low temperature. 



5. Asymptotic entanglement and twin oscillators correlations 

In order to have a general view of the persistence of both entanglement and twin 
oscillators correlations for a common bath, we consider the respective asymptotic 
states. The asymptotic entanglement for identical decoupled oscillators dissipating 
in a common bath was given in Ref. [16] as a function of the temperature and the 
squeezing of the initial state. We show the phase diagram for presence and absence 
of asymptotic entanglement in Fig. 4, in the weak coupling limit. 

We then estimate d in the asymptotic limit (t — > oo), assuming thermalization 
for the mode a;+ = {xi+X2)/v^ and free dynamics for the mode x- = (xi — a;2)/v2 
depending on the initial state. By explicitly writing down d and by replacing all 
the entries of the two-oscillator covariance matrix by their asymptotic value, we 
conclude that the variance of the difference between the occupation numbers is 




400 




100 200 300 400 



Fig. 2. Dynamics of logarithmic negativity (gray), quantum discord (black), and max[0, — (i/4] 
(light gray) (see the main text for details) for uji = L02 and a common bath. Left panel: T = 0.1a;i; 
right panel: T = cji. 
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500 



Fig. 3. Dynamics of logarithmic negativity (gray), quantum discord (black), and max[0, — (i/4] 
(light gray) for a common bath. Here, uJ2 = i^i, T = O.lwi and A : 



■.0.2uj'f. 



negative when the following inequality is satisfied: 



{xi){2{x', 



I) 



(P-)(2(P^ 



1) 



1 < (xi 



(/. 



(12) 



Since the observables of the non dissipating mode x_ are oscillating in time, in 
order to find negative values for d, we will take the minimum in a period. It fol- 
lows then, as for entanglement, it is equivalent to start from a two-mode squeezed 
state or two (separable) squeezed states. Finally we obtain the phase diagram rep- 
resented in Fig. 4, where entanglement and negative d (or negative P distribution) 
are obtained below the corresponding lines. We see that twin oscillators correla- 
tions are achieved only for very low temperatures, in contrast with entanglement 
for which the detrimental effect of temperature can be compensated by stronger 
initial squeezing. As a matter of fact, for r > 1 there are no twin oscillators unless 
the temperature is very low (T < 0.25aji). This critical temperature corresponds to 
2(4) = 1. 

An important remark is that the threshold curve for d < is not continuous. 
In fact, for r = we have d — for any T. On the other hand, for r — > O"*", the 
presence of twin correlations is determined by the sign of (x^) — (p^), and they 
only appear below a finite critical temperature. 

From these results, we learn that the presence of a decoherence-free channel 
allows the preservation of the quantum character of the state, which manifests 
itself both through the presence of correlations (entanglement and discord) and 
the negativity of d. This latter characteristic, signature of a negativity of the P 
distribution, is more fragile than the other quantum markers considered here. 

Entanglement and twin correlations have been analyzed also in Ref. [501 in a 
different case, for a pair of optical field modes obtained from parametric down- 
conversion when the input light is in a thermal and separable state. Even if the 
states obtained from (different processes of) thermalization of a squeezed state, 
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Fig. 4. Phase diagram for entanglement (logarithmic negativity) and twin correlations {d < 0) 
as a function of the temperature {T/uj-i) and initial squeezing, in the weak coupling limit (7 = 
2 X 10~^ /iTLOi ) and for identical (oji = tJ2) and decoupled (A = 0) oscillators. The asymptotic 
state is entangled or displays twin oscillators correlations in the lower areas, as limited by black 
and gray lines respectively. 

as considered here, and the thermal state after unitary action of parametric down 
conversion, as in Ref. [301 are different, we see that in both cases entanglement is 
more likely to be found that twin oscillators correlations. 



6. Conclusions 

Wc have studied the effects of dissipation on the quantum character of the state of 
two coupled harmonic oscillators. We discussed the system quantumness considering 
the dynamics of different indicators, such as entanglement, quantum discord, and 
twin oscillators correlations, a signature of negative values of the P distribution. 
We find that in general all these quantumness indicators vanish in the asymptotic 
limit of large times when they reach the equilibrium (apart from an exponentially 
small value of discord) unless a decoherencc-frce channel exists. In our model this 
channel is represented by the coupling of the two resonant (identical) oscillators 
to a common environment. Still, out of resonance between the two oscillators, it 
would possible, in line of principle, to rest ore t he noiseless channel by unbalancing 
the coupling of the oscillators to the bath.^^ We have also found that, whenever 
entanglement cannot be asymptotically preserved, its death time becomes similar 
to that of the variance of the difference between the occupation numbers (d) . 

After showing the effect of detuning between oscillators, (common or separate) 
environment, oscillators coupling and temperature through few examples, we ob- 
tained the phase diagram for asymptotic entanglement and twin oscillators correla- 
tions. In the case of identical decoupled oscillators, we have found that in order to 
find a negative value of d at very long time, the temperature needs to be smaller than 
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the one allowing for asymptotic conservation of entanglement. Moreover, the effect 
of temperature on entanglement can be compensated by stronger initial squeezing, 
being this not the case for twin oscillators correlations. So, this quantifier (d < 0) is 
even more fragile under dissipation than the entanglement itself. This observation 
could be more general as it is similar to what found in a rather different system in 
Ref.EOl 

This asymptotic phase diagram has been obtained for identical oscillators and 
a common thermal bath. Once the frozen degree of freedom disappears (either by 
detuning the oscillators or by introducing separate baths), the system is expected 
to asymptotically thermalize in all degree of freedom. Its quantum character is 
then generally lost, unless oth er rn echanisms are introduced, as for instance driving 
the system out of equilibrium.'^ As a matter of fact, twin oscillators correlations 
are mostly related to the presence and the direction of squeezing of the damped 
eigenmode (that is the sign of (a;^) — {p\))- 
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